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Group velocity and causality in standard resistive RMHD 
I. INTRODUCTION 

Recently, relativistic ideal MHD (ideal RMHD) equations have been used more frequently 
to perform numerical simulations of plasmas around black holesA^. The ideal RMHD nu- 
merical simulations showed that around a rapidly rotating black hole, anti-parallel magnetic 
field configuration, where the magnetic reconnection will be caused spontaneously, is formed 
naturally. Furthermore, they showed that much electromagnetic energy is transported into 
the region between the accretion disk and the black hole ergosphere, stored in the region, and 
then will be released due to the magnetic instability or magnetic reconnection. To perform 
calculations of the magnetic reconnection around black holes, we have to use the "resistive" 
RMHD equations, including resistivity, but no one succeeded in such a calculation yet. This 
is because we worry that "standard" equations of the resistive RMHD with the "standard" 
relativistic Ohm's law may have the causality problem, and "generalized" RMHD equations 6 , 
where the group velocity of the electromagnetic wave in the plasma with the plasma pa- 
rameter larger than unity is smaller than the light speed c, are very complex to perform 
numerical calculations. That is, with respect to the standard equations, we had a concern 
that the group velocity v g of the electromagnetic wave becomes larger than the light speed 
c and it destroys causality if we assume the group velocity means the propagation speed of 
transportation of information or energy-. For example, it is suggested that in a coordinates 
frame where the plasma flows with the velocity larger than the critical velocity c 2 /v g , the 
amplitude of the wave packet increases spontaneously-. This kind of instability contradicts 
to causality and the second law of the thermodynamics. On the other hand, using a linear 
theory of electromagnetic wave propagation, we can prove that no signal propagates with 
a velocity greater than the light speed in a plasma and such acausal phenomenon never 
happens as discussed by Jackson 3 (see Appendix N} . 

To solve confusion with respect to causality of the standard resistive RMHD equations, 
we calculated expected situations of such acausal phenomena (instability) using the standard 
resistive RMHD equations in a flat space-time. In conclusion, the numerical calculations 
showed that such acausal phenomena never happen. That is, in such cases that acausal 
phenomena are expected, the edge of the wave packet becomes dominant compared to the 
bulk part of the wave packet and the edge of the wave propagates within the light speed. 
This means that superluminal group velocity does not mean breakdown of causality as 
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discussed before^. That is, the standard resistive RMHD equations permit the superluminal 
propagation of the infinite wave packet but not superluminal information propagation, thus 
do not break causality. We can use the standard resistive RMHD equations within causality. 
It is noted that in a real plasma, group velocity of the electromagnetic wave is equal to or 
less than the light speed due to the inertia effect and momentum of electrons or electric 
current, which are neglected in the standard resistive RMHD. To treat these effects, we have 
to use "generalized" RMHD equations^. 

In this paper, we report numerical calculations of electromagnetic wave packets using 
standard resistive RMHD equations. In Section HIl we introduce the standard resistive 
RMHD equations and explain the causality problem suggested before. In Section IHIl we 
show a numerical simulation of the electromagnetic wave in uniform, rest, resistive plasma, 
and in Section IIVI we try to calculate the propagation of the electromagnetic wave packet 
in the situation where acausal phenomena are expected. In Section [V] with respect to one 
of the causality problems, we clarify we should use the head velocity as the velocity limit 
of the electromagnetic wave in plasma, not the group velocity based on the linear theory 
of propagation of electromagnetic plane wave. In Section Wl\ we discuss validity of the 
standard resistive RMHD equations for astrophysical relativistic plasmas. In Section IVHt 
we summarize the results of the numerical calculations and present some remarks. 



II. STANDARD RESISTIVE RMHD AND SUGGESTED CAUSALITY 
PROBLEMS 

The standard resistive RMHD equations with 4-form in Minkowski spacetime are as 
follows: 

W) = o, (i) 

d v T» v = 0, (2) 

d v *F" V = 0, (3) 

d v F» v = J", (4) 

U V F^ = = ri(J»-p' c U»), (5) 

where d u , U u , and J u are the 4-derivative in Minkowski spacetime, 4- velocity, and 4-current 
density, F^ v is the field strength tensor, *F^ U is the dual tensor of F^, and T^ u is the 
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energy-momentum tensor, 

T „u _ pg ^ + hUI *ijv + F» a F va - ^ u F kX F kX , (6) 

(see Jackson-, Weinberg 7 ). The scalar variables, density p, pressure p, enthalpy density h, 
charge density p' e = —U U J U are observed by the rest frame of the plasma. Equation ([5]) 
presents the standard relativistic Ohm's law. Through this paper, we use units so that the 
light speed is unity, c = 1. 

The 3+1 formalism of the standard resistive RMHD equations is written as, 

&1P \^ d 



dt ^ dx 



ir = -Es^ (8) 

3 

i = -E^e--w) (9) 

i 

f = -E^- do) 

E^. = 0. 

* = £^> ( 13 ) 



^ + £ = -( Jl - Pe^)> ( 14 ) 

where P l = T° l is the 3-momentum density, e = T 00 is the energy density, B l = | JE fc e^ k Fjk 
is the magnetic field, i? 1 = _F 0i is the electric field, 7 = U° is the Lorentz factor, v l = U l /j 
is the 3- velocity, and p e = J° is the charge density. Using linearized equations of Eqs. (fTOj) - 
( I14p . we have the dispersion relation of the electromagnetic wave propagating in a uniform, 
rest plasma: 

V 
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where u is the angular frequency and k is the wave number of the electromagnetic wave. 
Here, we assumed the resistivity is constant over the whole plasma. Then, the group velocity 
of the electromagnetic wave propagating in the uniform plasma is given by 

2 "* =>1. (16) 



We found that the group velocity is always larger than the light speed. If we assume that 
the group velocity is the propagation velocity of information, we can send an information 
faster than the light speed. That is, this information propagates outside of the causality 
corn t = ±x (Fig. Ufa)). Furthermore, when we consider a new inertia coordinates frame 
O' — x'y'z', which moves with the velocity Vq > l/v g observed by the plasma rest frame 
O — xyz, the signal seems to propagate toward the past (£' < 0, but not "absolute" past) 
(Fig. Etb)). If we combine this type of information transport in the two streaming plasmas 
with the relativistic velocity v > l/v g , we may be able to send the information to the 
absolute past. This is strange because if it is true, we can know the today's information 
(for example, about stock rate) yesterday. This is one type of causality problems discussed 
before. The other type of causality problem is about spontaneous decrease in the entropy. 
Fig. [D(a) shows the electromagnetic wave packet propagation with damping. The damping 
rate is approximately given by Qui = \/(2rj) from Eq. (|T5l) . when k > \/(2rj). Here, let 
me note that the electromagnetic wave and then wave packet do not propagate in the ideal 
MHD case (77 = 0) because the damping rate is infinite (Qu — > 00). This damping is 
due to the convergence of electromagnetic energy to thermal energy of plasma. During the 
damping of the wave packet, the entropy of the system increases. On the other hand, when 
we see the same wave in the new coordinates frame O' — x'y'z' (Fig. [H(b)), we find the 
electromagnetic wave propagates toward the left direction and grows up. The increased 
electromagnetic energy of the grown wave packet comes from the thermal energy of the 
plasma. This process shows the spontaneous decrease in the entropy. This is also a strange 
phenomenon. If it is true, using this effect, we can solve the energy problem of the world. 
As shown above, the causality problem is categorized into the following two types: 

1. Transport of information toward the absolute past. 

2. Spontaneous decrease in the entropy. 
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In this paper, we investigate these problems of the standard resistive RMHD using numerical 
simulations. Basically, we use the method developed by Watanabe et al.~ 



III. INFORMATION TRANSPORTATION 

A. Propagation of a simple wave packet 

We consider a simple wave packet propagation with group velocity v g in a rest, uniform 
plasma. Here, we take the parameters of the wave packet as, 

a = 50, k = 2vr, t) = oo = -^=(±1 - i) = ivre^ 2 , 

where a is the characteristic length of the wave packet, fco is the typical wave number, and 
Uq is the angular frequency. The initial conditions of the wave packet are given by 



p 


= Po 


= 1, 


V 


= Po - 


= 0.9, 


V 


= 0, 




E x 


= E Z 


= E X = By = 0, 


Ey 


= K 


= E cos27rxe' x2/2a2 , 


B x 


= By 


= B° x = By = 0, 


B z 


= B° Z 


= £ cos (2nx + -) e~ x2 l 2 ° 2 


J X 


= Jy 


— J° — J° — 


Jz 


~ J 


= 2 3 / 2 nE cos2nxe- x2 ^ 2 . 



(17) 



Through this paper, we set o = 50 and Eq = 0.01. We show the results of the wave packet 
propagation in Fig. [2J We clearly found the shift of the wave packet from the initial position 
with Aa; = 4.2 and Ax = 7 at t = 3 and t — 5, respectively. This means that the infinite 
simple wave packet propagates with the group velocity v g = 1.4. This confirms the results 
of Koide-. It is also shown that the wave packet damps very quickly. (Note that the scales 
of the vertical axes of the panels in Fig. [2] are quite different.) 
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B. Propagation of a wave packet with a head edge 

We consider a more realistic situation of information transport. In this situation, we 
can not use the infinite length wave packet at least without the head of the wave packet, 
because it takes infinite time to produce such a wave packet. Figure [3] shows the wave packet 
propagation which has the edge at the head x = 32, t = 0. The initial condition is given by 
multiplying the factor 



to each initial perturbation variable of the electromagnetic wave packet expressed by Eq. 



Fig. [3] clearly shows the bulk of the wave packet propagates with the group velocity v g = 
1.4, as shown in the case without the edge of the wave packet (compare the peak of the 
wave packet part at x = and x = 4.2 at t — and t = 3.0, respectively). On the 
other hand, at the front of the wave packet, a relatively smoothed new structure shaped 
like combined upside-down solitons (we name it "head structure", hereafter) appears and 
becomes dominant compared with the primary bulk part of the wave packet, because the 
head structure damps more slowly than the bulk of the packet does with a short characteristic 
wave length \ = 2ir/k = 1. The front of the head structure propagates with (less than but 
almost) the light speed. It should be emphasized that the head structure becomes dominant 
before the bulk part of the wave packet reaches to the initial position of the wave packet front 
edge. It is noted that even if the main part of the wave packet has an information, we must 
begin with creating the front of the wave, and then the velocity of the information transport 
never exceeds the head velocity of the wave, i.e., the light speed. It is clearly shown that 
the wave packet can't transport the information with the group velocity (which is larger 
than the light speed in the present case), and then the group velocity is not regarded as the 
velocity of the information propagation. As shown in Section |Vl the velocity limit of the 
information propagation is given by a variable called the head velocity. Now, the first type 
of the causality problem of the standard resistive RMHD equations, about the information 
transportation to the absolute past, is solved. 




(18) 



(HU) as 



Ei(x) = f{x)E°{x), B^x) = f(x)B?(x), and J^x) = f(x)J?(x). 



(19) 
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IV. SPONTANEOUS DECREASE IN ENTROPY OF WAVE PACKET 

We consider the second type of the causality problem in the standard resistive RMHD, 
that is, the suggested strange process of spontaneous decrease in the entropy of the electro- 
magnetic wave packet. Before we calculate the wave packet propagation in the new coordi- 
nates frame O' — x'y'z', where the electromagnetic energy of the wave packet is suggested to 
increase spontaneously, we have to consider the initial condition of the wave packet. When 
we use the approximation of the very narrow wave-number spectrum of the wave packet^, 
the perturbation variables of the wave packet, for example, the transverse component of the 
electric field, E y is given by 



EJx,t) oc — = exp 

yK ' ^1 + iDt/a 2 



1 

2a 2 1 + iDt/o 2 



e ~ik x-iuj t ^2Q) 



where uq = u(k ), v g = (du/dk)k=k , and D = (d 2 u/dk 2 )k=k "- As shown in Appendix 
|A| the profile of the initial condition (Eq. flA8j) ) in the new coordinates frame O' — x'y'z' 
monotonically decreases with respect to the coordinate x' . This is because the point at 
x' < 0, t' = corresponds to the past of the original coordinates frame, t < 0. In the past, 
the amplitude of the wave packet is (much) larger than the initial profile in the plasma rest 
frame. This electromagnetic wave with the monotonically decreasing profile in the frame 
O' — x'y'z' cannot be used as the transporter of information. This is because, in the limit of 
x' — > —oc, the amplitude of E' observed in the frame O' — x'y'z' is approximately given by 
\E' y \ ~ exp[— ('-/qVqD 2 / (r]a 2 ))x'}, from Eq. ( 1A14I) . which is significantly large over an infinite 
region. Here, vq is the relative velocity of the frame O' — x'y'z' and the frame O — xyz 
and 70 = (1 — Vq)^ 1 ^ 2 (see Eqs. ( 1A5I) and ( 1A6I) ). We cannot create such electromagnetic 
wave from a start source point of information in finite time and with finite energy. Then, 
to perform a simulation of the propagation of the signal with the wave packet, we have to 
introduce an edge at the rear of the wave packet. 

To avoid the Lorentz transformation at —x' ^> 1 of the new coordinates frame, we shift 
the wave packet with Ax = 32 toward the left side at t — (Fig. H]). We also confirmed the 
primary wave packet propagates with the group velocity v g = 1.4. We set the edge of the 
wave packet at the region —6 < x < 0, which is the relaxation zone of the step function at 
the wave packet edge. 

The initial condition is given by multiplying the following factor to each initial variable 
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of the electromagnetic wave packet as shown in Eq. ([191) ; 

[ (x< -I) 

f{x) = < 1 - (f ) 2 e (* 2 -' W (-/ < < 0) (21) 
1 (x > 0) 

where I = 6 through this paper. As the same of the front edge of the wave packet edge, 
the relatively smoothed structure appears around the wave packet rear edge and becomes 
dominant compared to the bulk wave packet (Fig. [5](a)). Fig. [5] (b) zooms up the edge of 
the wave packet and shows that the left edge propagates with the light speed. In this case, 
we can set the initial condition of the electromagnetic wave as shown in the first panel of 
Fig. El Because the left of the x'-axis of the new coordinates frame corresponds to the future 
of the primary coordinates frame O' — x'y'z', the bulk part of the wave becomes too small 
to see. Only the rear edge part of the wave can be seen. The edge part propagates toward 
the left direction with the light speed and damps (Fig. [6]). This is reasonable because the 
edge part propagates within the light speed in the plasma rest frame (Fig. [5]). Then, the 
second type of expected strange phenomenon of the causality break-down never happens. 



V. WAVE HEAD VELOCITY AS VELOCITY LIMIT OF INFORMATION 
PROPAGATION 

In addition to numerical results with the electromagnetic wave packet, we theoretically 
confirm that the group velocity does not give the limit of the propagation velocities of in- 
formation and energy in this section. With respect to the information velocity, we have 
confirmed numerically in Section UTTl that the group velocity presents the propagation veloc- 
ity of the wave packet without front or rear edge whose Fourier component is a Gaussian 
function with a narrow spectrum. While most part of the initial profile of the wave packet 
is concentrated in the region with the length of several a, the amplitude of the profile is 
finite over the whole range, — oo < x < oo. In the inertia frame, whose relative velocity to 
the plasma rest frame exceeds the inverse of the group velocity (v > l/v g ), the wave packet 
does not transport an information because its shape becomes a monotonic function and has 
no peak (Appendix . Then, we recognize such a wave packet which spread over the whole 
region — oo < x < oo is not accepted as an object of information transporter in relativity. 
Concerning the energy propagation velocity, in the first place, we note that we cannot define 
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the velocity of energy propagation because energy is a temporal component of 4-momentum 
and we cannot consider the transportation of energy only. In relativity, we have to consider 
transportation of 4-momentum, not only the energy component. In fact, the group velocity 
Vg = du/dk is not a scalar, much less 4- vector, and then the superluminal group velocity 
v g > c is not an invariant statement in any inertial coordinates frame. This contradicts to 
relativity. Then, we understand the propagation velocity of energy has no physical meaning 
in relativity. 

Instead of the propagation velocity of the wave packet without the edge (i.e., the group 
velocity), we can introduce the "head velocity", 



as the velocity limit of the propagation of non-zero fields. The head velocity t>h is a scalar if 
it is the light speed in a certain coordinates frame, and it is less than the light speed in any 
inertia frame if it is in a certain frame (Appendix [B]) . It has been already shown that this 
head velocity represents the information velocity limit in the dielectrics (see section 7.11 of 
Jackson-, and references therein). 

Here, we prove that the head velocity of the wave provides the velocity limit of information 
propagation also in a plasma. This proof confirms the results of Section II III with respect to 
the first problem of causality shown in Section [Til To investigate the propagation speed of 
the electromagnetic plane wave front, we consider an injection of the electromagnetic plane 
wave normally incident from the vacuum on the semi-infinite uniform plasma filling in the 
region x > 0, whose boundary locates at x = at t = (Fig. [7]). We assume the plasma has 
uniform resistivity rj, and the dispersion relation of the electromagnetic wave is presented 
by Eq. ([15]) . The index of refraction n(u) is given by 



According to Eqs. (7.124) and (7.125) of Jackson 3 , the amplitude of the electric field of the 
plane wave for x > is given by 



v h = lim — 



(22) 




(23) 




(24) 



where 




(25) 
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is the Fourier transform of the real incident electric field Ui(x,t). Suppose that the incident 
wave has a well-defined front edge that reaches x = at t = 0. Then, the incident electric 
field Ui(x,t) evaluated just outside the plasma at x = 0~ has the condition, 

Wi(0,t) = (t<0). (26) 

This condition, together with certain mathematical requirements that are "physically rea- 
sonable", is both necessary and sufficient to assure that A(oj) is analytic in the upper half 
of the complex uj plane. In this situation, the "physically reasonable" requirement is that 
«i(0, t) — > as t — > oo. This is true at least for the wave packet, while for semi-infinite 
continuous wave train at < t < oo, the physical reasonable condition is not satisfied*. 



For 



an example, we consider a simple case of 



Eoe-^coswot (f>0) 
Mi(0,t) = { , (27) 

(t < 0) 

where Eq, u, and ujq are the positive variables (when v — 0, it is not the case of "physically 
reasonable"). Then, we have 

, , . iEn uj + iv 

A(co) = — % 2 2. 28 

It has singular points at u — ±uo — ii> and a zero point at u — —iv (Fig. It clearly shows 
that there is no singular point in the upper half of the u plane. In the limit of \u\ — > oo, 
we have A — >■ then \A\ < ^3j- This inequality can be generalized for an arbitrary 
incident dumping wave train like a wave packet, if we take the constant E large enough. 

We formally evaluate the amplitude presented by Eq. f f24|) by contour integration. For a 
plasma case, the path of the contour integral includes the branch point of k(u) = u(u + -) 
at uj = 0. To avoid the multi-variable definition of k(u) around u = 0, we have to modify 
the integral path near the branch point u = so as not to cross the branch line between 
u = and u = — ^ as shown in Fig [91 Here, the radius of the half circle C e around uj = is 
infinitesimally small (e — > 0). With the definition of the head velocity (Eq. (122"]) ). the 
argument of the exponent in Eq. T2M becomes iuj(x — v^t) for large u, and the closing of 
the contour can be done in the upper half plane for x > v-^t, and the lower half plane for 



1 It means that the initial condition of the rear part of the wave very far form the front edge determines 
the propagation of the wave front. This case is out of our scope. 
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x < v h t. With respect to the case of x > v^t, this is because we have the inequality with 
respect to the integration along the large upper half circle Cr : z = Re 10 (0 < 9 < 7r), where 
R is the radius of the half circle Cr, 



c R 1 + n{u) 



A{uo)e ik{ul)x - iwt dLo 



M _ e —{x—Vtit)R\ 

£ Ea ( x -v>t)R ' (29) 



and then the integral along Cr vanishes in the limit of R — > oo. In Fig. [9j we can see the 
integration contour consists of C e , Cr, and I e< R for the case of x > v^t. With the analytic 
nature of n(u) and A(oj) in the upper half of the u Riemann plain, the whole integrand in 
the right-hand side of Eq. ( 12^1) is analytic there, and the Cauchy's theorem guarantees that 
the integral vanishes. We have thus proved that 

u(x,t) = (x>v h t) (30) 

provided that A(u) and n(u) are analytic for ^su > and n(u) — > 1/fh fo r 1^1 — ^ oo. 
Equation (130|) . together with Eq. (|26|) . establishes that no signal propagates with a velocity 
greater than the head velocity, in the plasma. This confirms the results of Section IIHI about 
the velocity limit of the electromagnetic wave packet, but not those of Section [IV] with 
respect to the causality problem of decrease in the entropy. In the plasma, the head velocity 
becomes the light velocity, — 



VI. USAGE OF STANDARD RESISTIVE RMHD EQUATIONS 

Now, we understand the standard resistive RMHD equations have no causality prob- 
lem, although the group velocity of the electromagnetic wave calculated by these equations 
exceeds the light speed. We showed that the head velocity of the wave determines the infor- 
mation propagation velocity limit as shown in the previous section. Then, we are free from 
the restraint of the superluminal group velocity in the standard resistive RMHD equations. 
Here, we consider the validity and limit of the standard resistive RMHD equations based on 
a more general point of view using the generalized RMHD equations. We already evaluated 
each term of the generalized general relativistic MHD (GRMHD) for the plasmas around 
several astrophysical black holes 6 , where we compare these terms with the V x B term. We 
compare the terms beyond the standard resistive RMHD equations with the resistive terms 
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in generalized relativistic Ohm's law, first, 

ne at \ne ) ne 2ne 

/ ^ + / _ A»j\ xB _ v[J _ , {1 + Qm (31) 

\ ne J \ ne J 

in flat space-time, which comes from Eq. (63) of Koide^. Here, n is the particle number 

density n = p/m, H^jt is the momentum density of the electric current, K is the stress 

tensor of charge and current, m is the characteristic mass of a particle m = m + + m_, 

a is the normalized reduced mass, a = m + m_/(m + + m_) 2 , A/i is the normalized mass 

difference Au = (m + — m_)/(m + + m_), Ap = p + — p_, and G is the thermal energy 

exchange rate from the negatively charged fluid and the positively fluid, where m± is the 

mass of positively/negatively charged particle, p± is the pressure of positively/negatively 

fluid in two- fluid model. The left-hand side of Eq. fl3Tl) presents the inertia of current, 

the first term of the right-hand side shows the transport of electric current momentum, the 

second term indicates the thermoelectromotive force, and the term np' c QU in the last term 

expresses the redistribution of the thermal energy due to the friction between the two fluids. 

We evaluate these terms beyond the standard relativistic Ohm's law except for the last 

term about thermal energy redistribution of friction. 

• Inertia of electric current: the left-hand side of Eq. ( T3T1) . 

(inertia of current) 1 1 d f\ih \ 1 u)/j,hJ u> 
(resistive term) i]Jnedt\ne J nJ {ne) 2 ^V e * ' 

where u is a characteristic frequency of phenomena, J is a characteristic value of 
the current density, and v e * is the Coulomb collision frequency between electron and 
positively charged particle (*=i: ion, *=e: positron). Here, we used the relations, 
rj = mv e */{ne 2 ) =>■ hv e * j '{ne) 2 , co p = ne/y^jih, and co c = enB/h. The Coulomb 
collision frequencies of electron-positron (z/ ee ) and electron-ion (z/ ei ) are evaluated by 

^LexlO-d)"' 2 ^)^], (33) 

„ d = 6.3xlO»(D" /2 (^)[s-], (34) 

where unit of T e /e is eV and unit of n e is m~ 3 (Appendix in Miyamoto-). Then, these 
collision frequency are comparable, z/ ee ~ 2.5f e j. 
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Thermoelectromotive force: the second term of right-hand side of Eq. (I3TI) . 
(thermoelectromotive force) 1 1 



(resistive term) 77 J 2ne 



|V(A//p- Ap)\ 



1 1 1 , A A , Afip - Ap /iWp 



-lA/xp - Ap| ~ " r P , (35) 

7/ J 2ne L 2n iy e *oo c 

where L is a characteristic length of the phenomena. 

Hall effect: the third term of right-hand side of Eq. ()3T1) . 

(Hall term) _ 1 A„, J x fl| ^ J_ A^ J5 = ^ (3g) 



(resistive term) r/J ne r/J ne is, 

where B is a characteristic value of magnetic field. 

In the generalized RMHD equations, the stress tensor is given by 

T ij = p ^ + hU i W + ^^(U i J j + J i W) + -^-rJ i J 1 + (— + —\ 5 ij — BiBj—EiEj. (37) 

ne (ne) 2 \ 2 2 / 

The momentum stress tensor of electric current is shown by the third and forth terms of the 
right-hand side of Eq. (13~T|) . The terms are evaluated by comparison with the kinetic stress 
tensor hWU^ as 

(current momentum tensor) 2uAhU J J ri B 2louj c 

< 2fj,— — ~ 2^— = — — . (38) 



(hydrodynamic stress tensor) U 2 neh Une U L r/ne u 2 

The magnetic Reynolds number Sm is defined by 

UL (nefUL ^UL fiu 2 p L u 2 p 
a M — — — r — S — A 4 5 l dy J 

where U is a characteristic values of 4- velocity \U\. Table 1 shows the values of these 
ratios for plasmas around several types of astrophysical black holes with the mass Mbh- 
Here, we assume that the plasma is electron-ion plasmas. We use the wavelength of the 
fastest-growing magnetorotational instability (MRI), Lqs, which causes the current sheet 
in the plasma disk as the characteristic length of phenomena around the black holes L. It 
is noted that the wavelength of the fastest-growing mode gives the maximum thickness of 
the current. Because in the linear stage, the magnetic field lines are bent with the scale of 
the wavelength and the current sheet is formed with the scale. In the nonlinear stage, the 
current sheet is pinched until the magnetohydrodynamic equilibrium between the pressure 
gradient and Lorentz force is achieved. It was shown in numerical MHD calculations (for 
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examples,— 1 ^). The wavelength of MRI at maximum growth rate is given by Lcs ~ 4y 
where f2 is the angular velocity of the disk around the black hole, f2 ~ a/ GMbhA' 3 > = 
2GMbh and va is the Alfven velocity, va = \j B 2 / p (G is the gravitational constant). It 
shows that the magnetic Reynolds number is much larger than the unity and the resistivity 
is negligible for the scale of current sheet (see also page 1468, right column of Koide-). 
With respect to Ohm's laws, we can neglect the inertia of current in all cases. Hall effect 
is negligible for supermassive black hole cases, while not for gamma-ray bursts (GRBs) 
and black hole (BH) X-ray binaries (including stellar mass black holes). On the other 
hand, the thermoelectromotive force would become significant compared to the resistivity 
in all situations. Because of the huge values of magnetic Reynolds number in Table 1, the 
electromotive force due to resistivity is very small compared to the ideal MHD electromotive 
force and then Hall effect is the same. With respect to the momentum transportation of 
current, it is negligible for all cases. In conclusion, we can use standard resistive RMHD 
equations with additional terms of thermoelectromotive force (and Hall term for GRB and 
BH X-ray binary.) 

When we consider the pair plasma, the normalized reduced mass and normalized mass 
difference become fi = 1/4 and A/i = from fi = m e /mi ~ 1/1800 and A/i = in 
the electron-ion plasma, respectively. Hall effect disappears and other special effects of 
generalized RMHD may be greater than the ratios of normal plasma by m[/m c ~ 1800. 
The conclusion with respect to the validity of the standard resistive RMHD equations is not 
changed. 

VII. SUMMARY 

We investigated whether the standard resistive RMHD equations have the causality prob- 
lem or not using the direct numerical calculations with the equations. There were two kinds 
of the hypothetical causality problem: (i) information transportation to the absolute past 
and (ii) spontaneous decrease in entropy. The numerical calculations showed that such phe- 
nomena do not appear even in cases where such phenomena are expected to happen. In the 
calculations, we used the electromagnetic wave packet with the characteristic length a = 50, 
characteristic wave number ko = 2ti in the plasma with resistivity 77 = l/(2 3 / 2 7r). These 
results suggest that the "group velocity" does not have physical meaning, such as the prop- 
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agation velocity of information. In fact, with respect to the first causality problem (i), we 
confirmed that the head velocity of the wave represents the velocity limit of the information 
propagation and the head velocity is the light speed in the plasma. 

With respect to the second problem of causality (ii), the calculation results suggest that, 
in general, acausal phenomena with a wave packet never happen because the wave packet 
without an edge in the plasma rest frame never be transformed to the wave packet with a 
peak by the Lorentz transformation with the relative velocity vq > l/v g , as shown in Fig. 
[TU1 (Appendix IA"]) . and the wave packet with an edge induces the (pair-soliton-like) head 
structure at the edge of the wave packet independently (or weakly dependent) on the bulk 
wave packet. The former shows that the concept of a wave packet is not consistent with the 
frame-invariant concept of relativity. Consequently, we understand that we have no causality 
problem in standard resistive RMHD as suggested by discussion with the superluminal group 
velocity of the electromagnetic wave^. This conclusion gives the guarantee of an usage of the 
standard resistive RMHD. This problem was the (mentally) biggest obstacle for application 
of the standard resistive RMHD equations to the global astrophysical phenomena. We expect 
that numerical simulations with the standard resistive RMHD equations to be more popular 
for analysis of plasmas around black holes, where the magnetic reconnection may play an 
important role in the release of the energy stored in the plasma near the black holes. 

Incidentally, it is noted that in the plasma, whose plasma parameter is much larger than 
unity, the group velocity is less than the light speech. To include this effect correctly, we 
have to include the inertia effect of the electron or electric charge/current. That is, we 
have to use generalized RMHD equations to consider the high frequency phenomena, where 
we cannot neglect the electron or current inertia, while it is not necessary for the global, 
relatively slow phenomena around the astrophysical black hole. 
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Appendix A: Initial condition of wave packet in a case of hypothetical acausal 
effect 



We show the profile of the electromagnetic wave packet propagating along the x direction 
becomes monotonic when we observe it from the coordinates frame to the relative velocity 
Vq > 1/vg to the rest frame of the plasma, O — xyz. Here, in the frame O — xyz, we assume 
that the Fourier spectrum of the wave packet is Gaussian, 

o 



F(k) 



'2n 



-e 2 



(fc-fco) 2 



(Al) 



where a corresponds to the characteristic length of the wave packet and ko is the typical 
wave length of the wave packet. Furthermore, we assume a is large enough and the Fourier 
spectrum is significant only near k ~ k . We have the complex perturbation variables of the 
electromagnetic wave packet as 



E y (x, t) oc 



1 



y/l + iDt/a'- 



cxp 



2o 2 1 + %Dt/a 2 _ 



-ikQX—iujQt 



(A2) 



where ujq = uj(k ), v g = (du/dk)k=k , and D = (d 2 u/dk 2 )k=k —- The physical solution is 
given by the real part of the complex perturbation, E y . The amplitude around any point 
x = Xq at any time t = t is given by, 



E y (x, t) oc 



y\ + iDt/a 2 \ 
1 

VI + (Dt/a) 2 



exp 



cxp 



3? 



1 (x — v g t) 2 
2a 2 1 + iDt/a 2 
1 (x - v g f 



,(SJu)o)t 



2a 2 1 + (Dt/a) 



(9fu)o)t 



(A3) 
(A4) 



The Lorentz transform from the original coordinates frame O — xyz to the new coordinates 
frame O' — x'y'z' with the relative velocity v = v x, 



t' = 7o(* - v x),x' = 7o(x - v t), y' = y,z' = z, 



(A5) 



or 



t = 7o (t' + v Q x'),x = 7o(x' + v t'),y = y',z = z', (A6) 

where 70 = (1 — v 2 ,)^ 1 / 2 . Then, the amplitude profile of the initial electric field of the wave 
packet observed by the new coordinates frame O' — x'y'z' is 

v k 



\E'\ = \ l0 (E y -v B z )\= l0 



LJ 



\E„ 



1 



1 + {D l0 v x'/a 2 ) 2 Y^ 



exp 



1 — v g v ) 2 x' 2 



2a 2 1 + (D~/ v x'/a 2 ) 2 _ 



D 701)0(3^0)2;' 



(A7) 
(A8) 



17 



Group velocity and causality in standard resistive RMHD 



where we assumed v ^ u) /k . Here, we have the dispersion relation of the electromagnetic 
wave in the plasma with the uniform resistivity rj, given by Eq. ffT5]) . Then, we have 



CJ 



2rj k 



D 



\AW - 1 
1 ( 277 



2 V 2 \v^2W 



< 0. 



Here, we have the relation, 



D 

V 



1 /Vg 

2 V 77A; 



(A9) 
(A10) 

(All) 
(A12) 



We note that \E y \ becomes infinite when x' — > — oo because $tuo = —1/t] < 0, D ^ 0, and 
Wo / 0. This means that this wave packet keeps its finite profile only in the original frame 
O — xyz. In other coordinates frame O' — x'y'z' with the finite relative velocity to the rest 
frame of the plasma, O — xyz, the electromagnetic wave can not be recognized as a wave 
packet anymore. Then, we understand that the concept of "wave packet" is not invariant 
item with respect to the Lorentz transformation. 

To show the monotonicity of the initial perturbation in the new coordinates O' — x'y'z', 
we evaluate the following variable proportional to the logarithm of the profile of the pertur- 
bation, 



1 - v g v ) 2 a 2 



-x 



J2 



tX 



2rj 2[a* + (D l0 v x') 2 Y~ i 
Using a new variable £ = — 10 ^l D x', we have 



log{a A +( lo v Dx') 2 } 



(A13) 



V 



(A14) 



We evaluate the slop of L, 

dL _ D( l0 v ) 2 
«9£ 



+ {VgV 



2e 



2\2 



(A15) 



7i l a 2 i + ' '' (l+£ 

When £ is greater than zero, it is clear that dL/d^ is negative. When £ is negative, dL/d^ 
is negative if 

3 

„ „ „4n z 

< 0, 



An 2 

(v g v - 1) 2 (1 - v 2 ) + 2w > 2 - l) 3 -^ 



18 



Group velocity and causality in standard resistive RMHD 



where we use inequalities, |2f/(l + £ 2 )| < 1 and |2f/(l + £ 2 ) 2 | < (V3/2) 3 . When v > l/v g 
and v g > 1, we have 

(Vo-D 2 (l-^)<(«r-l) 3 . 

3 



and 



2/3 , 
vJ — 1 



< (2cos4tt/9) 



3/2 (1-2cos47t/9) 3 



1)3/ 2 v ' 7 (1-8cos 3 4tt/9) 3 / 2 ' 

where 2 cos4tt/9 = 0.3472 • • • . Then, dL/d^ is negative if 



1|<2, (, 2 -1)3 



< 2v {vl - If 



{vl - If I 2 



1 _ I ( yl 1 x 0.35 3 / 2 (l - 1/3) E 



!)3/2 
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4 \ 2 

4r/ 2 



4?7 2 



< 0. 



On the other hand, because we neglect the third order terms of | A; — A; 1 < l/cr with respect 
to uj, we already assumed that 



1 d' A uj 1 
3! dfc 3 <7 3 



2! 9fc 2 a 2 



1 <9L> 



a 



(A16) 



where D = d 2 u/dk 2 = -4r]/[(2r]k) 2 -1} 3/2 . When we use the relation ((2r)k) 2 -I)- 1 = v 2 -l, 



we have 



4r]k(v 2 — l)r) 



< 1 



Using the relation, ^ = 2(t> 2 — l) 1 / 2 , we have 



77 

- < 



a 27^2 _ l)V2 



and then we obtain 



< 



< 



Then, we find 



4vf 

a 2 " 2v 2 g {v 2 g - 1) " 7j3(l - 1/^2)3/2 - ( v 2 _ !)3/2' 

4r? 2 ° 



< 



(7j2 - 1)3/2 



x 0.989993. 



This confirms dL/dl; and then dL/dx' are always negative. This calculation shows that the 
slope of the transformed shape of the original wave packet into the new coordinates with 
vo > l/v g is monotonically decreasing function, and its shape is drastically different from 
that set in the plasma rest frame. That is, the concept of "wave packet" with infinite length 
is not invariant for the Lorentz transformation, and contradicts with relativity. 
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Appendix B: Lorentz invariant property of head velocity 



We show the frame transport property of the head velocity of a wave train. We consider 
two different inertial frames x M and where the metric is given by 



ds — jj^jjdjtic diQij — T~j pQ-d3j doc • 



(Bl) 



Then, we have the relation with respect to the transportation coefficient between the two 
frames, 
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The transform of the 4-wavenumber k^ = (u, k) is given by 



to 



dx ' 
dx° 

dx 1 ' 



UJ 



dx° 



UJ 



dx ' 
dx 1 

dxi 



k' 



The wave head velocity in the coordinates frame O' — x'y'z' is calculated 
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where we define n J = fc'/A;. Using Eq. (1B2j) . we derive 
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3 / \ 3 

Then, if < 1, we have < 1. The equality stands only if = 1 in a certain frame. 
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FIG. 2. Simulation of propagation of a simple wave packet with a narrow Fourier spectrum. 
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FIG. 3. Simulation of a wave packet with a head edge. The small boxes in the top and bottom 
panels show the zoom-up of the profiles of the head edge of the wave packet in the range 30 < x < 
34, -10" 4 < E y < 10" 4 and 35 < x < 39, -1(T 12 < E y < 10" 12 at t = and t = 5, respectively. 
The latter small panel shows that the head edge locates at x = 37 at t = 5. 
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FIG. 4. Simulation of propagation of wave packet with a very narrow Gaussian Fourier spectrum. 
The initial condition of this wave packet shifts toward the right direction by Ax = 32 from the 
calculation of Fig. (TJ 



23 



Group velocity and causality in standard resistive RMHD 



(a) 



10" 



-10" 



1 1 1 1 1 1 II 1 1 II 1 111 1 1 1 1 ' 

.III III,,, 



I ^w\AJ 

-10 



5X10" ! 



2X10 



-2X10" 
6X10" 



-8X10" 



10" 



10" 



10" 



10" 



10" 



10" 



12 -10 



-12 



-10 



-12 



-10 



-6 

X 



t=0 




-10 10 20 30 40 50 60 70 




-10 10 20 30 40 50 60 70 



t=5 



-10 10 20 30 40 50 60 70 

X 

(b) 




FIG. 5. 

Simulation of propagation of a wave packet with a back edge in the plasma rest frame O — xyz. 
(a) Whole profiles of the wave packet, (b) Profiles of the wave packet around the rear edge. The 
small box in the top panel shows the rear edge of the wave packet in the range — 10 < x < 0, 



-5 x 1(T 5 < E y < 5 x 1(T 5 . 
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FIG. 6. Simulation of propagation of a wave packet in the coordinates frame where the plasma 
flows toward the left direction with vo > l/v g . This coordinates frame corresponds to the new 
frame O' — x'y'z 1 . 
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FIG. 7. Injection of electromagnetic simple plane wave from the vacuum to the semi- infinite 
uniform plasma at x > 0. 
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FIG. 8. Singular and zero points of A{uj) on the u complex plane. 
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FIG. 9. Branch points and a branch line of k{u) on the complex u plane and the path of the 
contour integral for the amplitude of the electric field (Eq. ([24 
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FIG. 10. Corrected image of the panel (b) of Fig. [TJ In the frame O' — x'y'z', the profile of the 
electromagnetic wave becomes monotonic and transports no signal because it has no peak. 
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TABLE I. Microscopic variables (ut u% ',•••) and characteristic scales of phenomena of plasmas 
around black holes to examine validity of the standard resistive RMHD equations. 
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a Data are from McKinney 12 

b The minimum values of L and r are estimated by the thickness (Lqs) and the light transit time (Lcs/c) 
of the current sheet caused by magneto-rotational instability (MRI) in the disk. 
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